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ABSTRACT 



We report on a comprehensive set of observations of Gamma Ray Burst 050716, 
• detected by the Swift satellite and subsequently followed-up rapidly in X-ray, optical 

CU an d near infra-red wavebands. The prompt emission is typical of long-duration bursts, 

with two peaks in a time interval of T90 = 68 seconds (15 - 350 keV). The prompt 
£^ ■ emission continues at lower flux levels in the X-ray band, where several smaller flares 

^ I can be seen, on top of a decaying light curve that exhibits an apparent break around 

03 220 seconds post trigger. This temporal break is roughly coincident with a spectral 

break. The latter can be related to the extrapolated evolution of the break energy in 
the prompt 7-ray emission, and is possibly the manifestation of the peak flux break 
. frequency of the internal shock passing through the observing band. A possible 3 a 

' change in the X-ray absorption column is also seen during this time. The late-time 

afterglow behaviour is relatively standard, with an electron distribution power-law 
index of p = 2; there is no noticable temporal break out to at least 10 days. The 
broad-band optical/nIR to X-ray spectrum indicates a redshift of z > 2 for this burst, 
with a host-galaxy extinction value of Eb-v ~ 0.7 that prefers an SMC-like extinction 
curve. 
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1 INTRODUCTION 

One of the main goals of the Swift mission (|Gehrels et alj 
|2004| ) is to obtain early-time information for gamma-ray 
burst (GRB) afterglows in X-ray and optical/UV wave- 
lengths, using its rapid and autonomous slewing capability. 
In addition, longer term monitoring, especially in X-rays, 
has become more feasible than previously, with a dedicated 
facility like Swift. 

Results of the early observations include the rapid de- 
cline of (X-ray) after g lows in the first tens of minutes 
i|Tagliaferri et all 120051 : iGoad et al.l 120061 ). An explanation 
for this behaviour is prompt e mission seen from angles fur - 
ther away from our line of sight (|Kumar fc Panaitesculfe oOO). 



It is likely that emission from the afterglow itself has not 
risen enough so early to contribute significantly to the mea- 
sured flux, although there are some exceptions where the (X- 
ray) flux of the late-time exte rnal forward shock i s already 
visible from early times on feg lO'Brien et~aLll2006l ). A thor- 
ough understanding of this phenomenon requires detailed 
descriptions of the light curve (and broad-band spectral) be- 
haviour early on, preferentially with as little as possible in- 
terference from other phenomena such as X-r ay flares. These 
X-ray flares are seen in a number of bursts feg lBurrows et all 
l2005al ). and are now general ly interpreted as continuing ac- 
tivity of the in ner engine l|King et alj 120051 : but see also 
iPiro et al. I l2005l , who interpreted the X-ray rebrightenings 
as the result of the beginning of the afterglow). Late time 
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monitoring allows the construction of multiple epochs of 
combined near infra-red/optical/UV and X-ray broad-band 
spectra (occasionally including radio data), building a more 
complete picture of GRB afterglows than before. 

We present here a full analysis of GRB 050716. This 
includes the prompt (burst) emission in gamma-rays and 
early X-ray data, for which we performed time-resolved spec- 
troscopy. We present the late time X-ray behaviour of the 
afterglow, as well as late time optical observations, which 
we then combine to form a broad-band spectrum. The ob- 
servations and their analysis are described in Section [2] and 
in Section [3] we discuss the results of the observations. In 
Section [4] we summarise our findings and draw our conclu- 
sions. 

In the following, 1 a errors are used except where noted. 
The temporal and spectral power law indices a and f3 are 
defined by F oc t -c V -/3 , and the photon index r = 1 + /3. 
All Swift data have been reduced using the Swift software 
version 2.4 within the HEAsoft software (version 6.0.5), and 
the corresponding CALDB files. 



2 OBSERVATIONS AND ANALYSIS 

GRB 050716 was detected by the S wift Burst Alert Tele- 
scope (BAT, iBarthelmv et all l2005h on July 16, 2005, at 
12:36:04 UT. All times mentioned in this text are relative 
to this To, except where noted. The spacecraft slewed and 
started observing in the X-ray and optical /UV bands 96 
and 99 seconds later, respectively. The XRT dBurrows et al.l 
2005b) was able to immediately locate a position for the X- 
ray afterglow on-board, allowing ground-based telescopes to 
perform rapid, deep follow-up observations. 



2.1 BAT analysis 

BAT data were reduced using the BAT software tools within 
the Swift software package. A PHA to PI energy conver- 
sion was performed first with the bateconvert tool, which 
ensures the conversion is quadratic, rather than linear. We 
subsequently corrected the data for hot pixels. An additional 
systematic error correction was applied before spectral fit- 
ting. The 15-350 keV light curve is shown in Figure [1] 

The light curve shows an increase similar to the typ- 
ical Fast-Rise Exponential-Decay (FRED) behaviour seen 
for many other bursts, with a peak around 9 seconds af- 
ter trigger and the start of the rise of the light curve some 
30 seconds earlier. A notable difference is that the rise is 
not as steep as usually seen for a FRED. A similar FRED 
light curve is superposed on the tail of the first one, pro- 
ducing a second peak 39 seconds after the trigger. After 80 
seconds, most of the prompt emission has disappeared in the 
15-350 keV energy band. Tgo in this energy band for this 
burst is 66 ± 1 seconds. A 4-second binned light curve shows 
some emission still apparent past 80 seconds, up to 200 sec- 
onds, which is mostly found in the lower energy bands. 

We have subdivided the gamma-ray data into 4 tempo- 
ral sections and obtained spectral fits for each of these. The 
spectral shape was fitted with a power law, which results in 
an acceptable fit. A power law with an exponential cut-off, 
however, provides a better fit, while a fit with a Band model 




Time (seconds since trigger) 

Figure 1. Top panel The BAT 15-350 keV mask- weighted 
light curve, with the two peaks around 9 and 39 seconds. Points 
are 0.5-second binned light curve, while the overplotted line is the 
4-second binned light curve. The latter shows some evidence for 
possible emission after 80 seconds. Bottom panel The spectral 
evolution of the prompt emission in the 15—150 keV range; shown 
is the photon index for a power law model with an exponential 
cutoff. 

ijBand et al.lll993l ) does not improve the fit compared to ei- 
ther one of those models. The BAT spectral energy range 
and low number of energy bins do not allow an accurate de- 
termination of the cut-off energy E c in the cut-off power law 
model (for the first section, it is in fact a lower limit), but 
there is an indication that the spectrum softens, with E c 
gradually becoming lower, and the photon index V chang- 
ing from about 0.4 to 0.9 during the prompt gamma-ray 
emission, over a period of ~ 50 seconds. 

To verify this spectral evolution, we have fitted the later 
time sections (C & D) with the model parameters (other 
than the normalisation) fixed at those obtained from the 
earlier sections (A & B). By then freeing the parameters and 
applying an F-test, we find that there is indeed a change in 
the spectrum at a 99.9% level. Performing a similar test to 
the first and last sections (A & D), with the results obtained 
from the fit to the total section (-20 - 80 s), also indicates 
a changing spectrum at a 99.9% confidence level. 

Such softening behaviour of the prompt sp ectrum has 
been commonly found before in other bursts (|Ford et alj 
1995). The results are detailed in Table[T] and include the re- 
sult of a spectral fit to the whole time span. Figure Q] shows 
the spectral evolution of the prompt emission between 15 
and 150 keV as well. 

2.2 XRT analysis 

XRT data were reduced from level 1 to level 2 with the 
Swift software task xrtpipeline. Data obtained in Win- 
dowed Timing (WT) mode and Photon Counting (PC) mode 
have been used in all our following analysis. Counts were 
grouped by 20 per bin for the spectra. WT mode data for 
the light curve have been grouped by 40 counts per bin, and 
PC mode light curve data have been grouped by 20 counts 
per bin. 
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Table 1. Results of spectral fits to four sections of the prompt 
gamma-ray emission, as well as that of the combined fit. 



Section Time T cutoff energy x 2 

(sec) (keV) 



A 


-20 


- 


0.42 


+0.21 
-0.60 


221 (> 6.7) 


1.05 


B 


- 


20 


0.74 


+0.27 
-0.29 


125 i^ 8 


0.86 


C 


20 - 


40 


0.64 


+0.40 
-0.45 


105 tlf 


0.86 


D 


40 - 


80 


0.87 


+0.33 
-0.36 


65.6 till 


0.95 


A - D 


-20 - 


80 


0.85 


+0.21 
-0.23 


150 +lf 


0.76 



We have subdivided our analysis into an early-time part 
and a late-time part. We designate the early-time part as 
the section of the X-ray afterglow that shows flares (proba- 
bly related to inner engine activity, see also the discussion 
below), while the late-time part shows a relatively smooth 
decay and is presumed to be the afterglow emission from 
a forward shock. The data for the early part extends up 
to almost Iks, with most of the data having been obtained 
in WT mode. The late-time part starts some 4ks after the 
trigger and extends out to almost three weeks. 

2.3 Early XRT data 

WT mode observations started at 103 seconds, and extended 
up to 516 seconds. After this, the count-rate was low enough 
for the instrument to automatically switch to PC mode. WT 
mode data were extracted using a rectangular region with a 
length of 93" centred on the source and along the readout 
direction, with an equally sized region off-source serving as 
the background determination. PC mode data have been 
extracted with a 47" circular aperture instead of the default 
71" because of a contaminating close-by source. The first 
few 100 seconds of PC mode data are piled-up and we used 
an annular extraction region, with a 12" inner radius and a 
47" outer radius. In addition, some bad columns are located 
near the source centre. We used default grades for WT mode 
(grades to 2), while we used only grade for the first few 
100 seconds of PC-mode observations. 

To determine the combined correction factor for the an- 
nular aperture and bad column loss, we modelled the PSF 
using ximage, and calculated the ratio of the integrated re- 
sponse between an unmodified PSF and that masked out 
by an inner 12" circle and the bad columns. In the same 
manner, we correct for the expected loss of using only a 47" 
circular aperture. The combined correction factor is 4.47. 

2.3.1 Spectral evolution 

Similar to the BAT observations, we find evidence for a 
change in the spectral slope over time. Modelling the spec- 
trum with a simple power law, we find that the X-ray spec- 
trum starts with V = 0.9 ±0.1 (0.3 - 10 keV), and the slope 
then evolves to T = 2.0 ±0.1, with a change between 250 
and 400 seconds post trigger. 

If a break is expected in the X-ray band, a broken 
power law model for the spectra around this time is pre- 
ferred. Fits with broken power law models, however, do not 
show a significant improvement in the fit, and we have there- 
fore used the results from the single power law fits above. 

Because there is still some flux in the BAT light curve 



coincident with the first 100 seconds or so of WT mode data, 
we have performed a joined fit to the BAT and XRT data 
between 100 and 200 seconds post trigger. A single power 
law fit has a Xred°f 1-98 (degrees of freedom (DOF) = 61), 
but a Band model does a much better job (x? c d = 0.90, DOF 
= 59), and results in Eq = 12.2 jl^" \ keV, or a peak energy of 
E p = 13.1±| | keV. The spectral indices are a = -0.93t ;?f 
and (3 = -2.56ti.Ti- 

Because of the low number of counts in the BAT data, 
the fit is dominated by the XRT data and the outcome 
should be taken somewhat cautiously. 

2.3.2 Light curve 

Since there is a change in the spectral slope halfway 
through the light curve, we cannot simply fit the count-rate 
light curve: a flux calibrated light curve has to be used. To 
derive the count-rate to flux conversion factors, we have used 
the spectra between 103 and 223 seconds and between 263 
and 513 seconds (see also Tabled, which provide the con- 
version values for the first and last part of the light curve, 
respectively. For the middle section we have used the av- 
erage of these values: the rapid change in spectrum does 
not allow for a good spectral fit. In all our following fits, 
we have modelled the two most pronounced flares with two 
simple Gaussians. We also left out the last data point of the 
light curve: it is likely that this point is the onset of the 
later-time light curve (see Section or possibly even the 
start of another flare. 

A fit to the flux light curve with a single power law re- 
sults in Xrcd = 2.06 (DOF = 123), while a broken power law 
fit gives x? cd = 1-83 (DOF = 121). Neither fit is good, but 
the broken power law does provide a significant improve- 
ment. A fit with a smoothly broken power law (either with 
the sharpness of the break fixed or free) hardly improves this 
(x?od = 1-78 for both, DOF = 121 and 120, respectively), 
and an exponential fit is worse (x? c d = 2.07, DOF = 123). 

The broken power law fit results in decay parameters 
qi = 0.9lt°;og and a 2 = 3.78± ;^. The temporal break is 
found at 220^| seconds, just before the time of the spectral 
break. This coincidence strengthens the suggestion that a 
broken power law is a reasonable model for the underlying 
light curve. See Section [3. II for a more detailed exploration 
of this break. 

The light curve and the spectral evolution are shown in 
Figure [5] As commonly seen in other early X- r ay light curves 
(eelBurrows et al.ll2005al : iNousek et al]|2006l : lO'Brien et ail 
2006), this one exhibits a few flares, which most likely in- 
dicates activity of the burst engine itself. This would agree 
with the possible emission seen in the BAT 15-350 keV range 
between 100 and 200 seconds. Also visible is the deviation 
of the underlying decay from a power law. 

2.3.3 X-ray absorption change? 

For the combined WT-mode data, we find no evidence for 
excess Nn above t hat of the estimated Ga lactic column 
of 1.1 x 10 cm" 2 jDickev fc Lockman|[l990l ). Subdividing 
the WT-mode into a pre- and a post- spectral break sec- 
tion indicates, however, a modest amount of excess absorp- 
tion at early times, while no excess absorption is measured 
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Time (seconds post TO) 

Figure 2. Top panel 0.3-10 kcV early-time light curve. The last 
three points are PC mode data, the previous data points consist 
of WT mode data. A broken power law fit with two Gaussians 
is shown as the continuous line, with its break time around 200 
seconds post trigger (the Gaussians account for the two main 
flares, to obtain a better overall fit; the broken power law without 
Gaussians is shown as the dashed line). The last data point has 
not been fit, since this likely belongs to the later-time light curve. 
Middle panel The fit residuals, expressed in sigma deviation. 
Even with the Gaussian fits to two flares, the overall light curve 
shows significant deviations from the fit. Bottom panel The 
evolution of the photon index Y of the X-ray spectrum. 



past the break: Nh = 2.0 (±0.2) ■ 10 21 crrT 2 and 7V H = 
1.1(±0.2) • 10 21 cm~ 2 for the total Nn at the two respective 
epochs. While the observed change in photon index could 
correlate with this change in Nn , a contour plot (Figure [3]) 
indicates it is likely that the column density did change. 
Something similar has been seen, for examp le, in the early 
(prom pt) X-ray emission for GR B 000528 (jFrontera et al.l 
|2004 ). or claimed for GRB 050730 (IStarling et al.ll2005l ). An 
F-test for a model with A*h fixed at the later time value 
(see also Section l2.4p , compared to a model with A'h free 
to fit, gives a probability of 6.7e-5 for the two models to be 
consistent. See the second part of Table [2] for more details. 



2.4 Late time XRT observations 

The XRT data after 4ks are not piled-up, and were extracted 
using a circular aperture of 47"radius and the default grades 
(0 to 12). The resulting 0.3-10 keV light curve follows a 
smooth power-law decline, with a power-law decay index of 
0.99. The full light curve is shown in Figure [4] and includes 
the early time light curve as well. There is no indication of a 
shallow part (a < 0.5) in the light curve between the steep 
early decay and the later tim e decay, as often se en around 
this time in the light curve (eg lNousek et al]|2006l) ; this may 
have been missed if it was relatively short and between the 
early and late-time part. There has been, however, no need 
to include such a section in our fits to obtain good results. 
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Figure 3. Contour plots for the 1-cr, 90% and 99% combined con- 
fidence intervals for the photon index T and Nn. The contours 
upper-left are for the WT-mode data after the observed spectral 
break (303-513 seconds), the contours in the lower-right corner 
are pre-break (103-253 seconds). The middle part (253-303 sec- 
onds) is left out, since here the change in spectral break is most 
pronounced, indicating T is highly variable. Note that T for the 
first data section is slightly higher than the previously quoted 0.9, 
which is for the very first part of WT-mode data only. 



There is also no indication of a late-time break before the 
light curve is lost in the background noise, up to 10 6 seconds. 
A spectral fit to the data between lOks and lOOks, with A^h 
fixed at its Galactic value, results in a good fit (Xicd = 0.54), 
with a photon index of T — 2.01^ 



I +0.11 

-o.io- 



2.5 Optical observations 

Optical observations were automatically performed with 
the Swift Ultra V iolet and Optical Telescope (UVOT, 
iRoming et al]|2005l ) in the UBVW1 M2 W2 filters and with- 
out a filter. The prompt dissemination of the position also 
triggered the Faulkes Telescope North (FTN), which began 
the first observations of its automatic follow-up procedure 
l|Guidorzi et al.ll2006l ) 243 seconds after the Swift trigger. 
Finally, near infra-red observations at the United Kingdom 
Infrared Telescope (UKIRT) started 56 minutes post trigger, 
in the JHK bands, and Gemini NIRI if-band observations 
were obtained 22 hours after To. The log of our observations 
can be found in Table [3] 

A faint, fading source was detected in the UKIRT data, 
just outside the 90% XRT error c ircle, and identified as the 
IR afterglow (jTanvir et al.ll2005l ). The IR data also gives 
the best position for the GRB and its afterglow. This is 
RA = 22:34:20.73, Dec = +38:41:03.6 (J2000), with an esti- 
mated error of 0.4". This position is derived from the NIRI 
if-band image, calibra ted with respect to the 2MASS survey 
l|Skrutskie et al.ll2006h . 

A single power law fit to the three A'-band points results 
in a power law decay index of a a = O.SOIq'Jy, with Xicd ~ 
0.689. While different from the X-ray decay slope at these 
times, the paucity of data points does not allow a strong 
quantification of this difference. A fit with a Q fixed at ax = 
0.99, however, results in Xvcd = 2.58. 

We calculated the JHK spectral index by shifting the 
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Table 2. Results of spectral fits for the XRT WT-mode data. The second half of the table details the results on variation of the column 
density. See also Figure[3] A simple power law was used, with the Teubingen-Boulder model l|Wilms et al,[ 2000, tbabs) for the absorption. 



time T N H x? od 

(sec) (10 21 cm~ 2 ) 



103 - 


123 


0.95 ±0.12 


-1 oy+0.73 
J-3' _o.66 


0.509 (12.2/24) 


123 - 


143 


i O1+0.15 
± - o± -0.14 


-, oy + 0.74 


0.899 (15.3/17) 


143 - 


163 


1.26 ±0.14 


i 1O+0.66 
1-1S-0.61 


0.518 (9.84/19) 


163 - 


183 


i 1O+0.13 
1 - i °-0.12 


i 11+O.6I 
J " 11 -0.S7 


0.698 (16.8/24) 


183 - 


243 


1.23 ±0.11 


i kk+0.42 
1 -" J -0.39 


0.801 (25.6/32) 


243 - 


363 


1.66 ±0.07 


i ci +0.25 
1 - O1 -0.23 


1.09 (63.3/58) 


363 - 


403 


2 lfi+ - 23 


1 8R+ ' 56 


1.12 (13.5/12) 


403 - 


513 


2 05+ 17 


1 nq+ ' 39 
i.ua_ 36 


1.09 (22.9/21) 



Nh variation 

103 - 223 1.26 ±0.05 2.04j^;?,J 0.890 (111.2/125) 

263 - 513 l-86to.o? L13 ±ai8 L053 (73.69/70) 




10'' 10" 10' 1 

time (seconds since burst) 

Figure 4. The complete X-ray light curve between 0.3 and 10 
keV. The continuous line is modelled using a broken power law 
and 2 Gaussians for the early time light curve (see Section 12.31 1 , 
and a single power law added to that for the late-time light curve. 
The latter has a decay index of 0.99 ± 0.02. Bottom panel The 
residuals of the fit, expressed in sigma deviation. 



data to a common epoch using the previously calculated 
decay in dex, and correcting for the estimated Galactic ex- 
tinction (|Schlegel et al.| [l998). Co nversion to flux was per- 
formed using the Veg a fluxes from Fukugita et alj l|l995t ) for 
the optical and from iTokunaga fc Vaccaf i 20051 ') for the nIR 
filters. The fit shows the spectrum to be inconsistent with 



a power-law spectrum (/3 = 2.9, but Xrcd = 4.2), indicative 
of intrinsic reddening. We address this in Section 



3 DISCUSSION 

Here we discuss the results obtained above. A short sum- 
mary of the results is given in Table [4] 



3.1 Early emission 

GRB 050716 shows characteristics as seen in previous 
(Swift) bu rsts, both in its prompt emission and the (early) 
afterglow l|Nousek et alj 120061 ; lo'Brien et"all l200d ). The 
flares, combined with the multiple peaks in the gamma- 
ray emission, imply sev eral shocks, likely indicative of pro- 
longed engine activity IjKing et al.1 120051 ) . The early X-ray 
data show a power-law decline, that appears to be unre- 
lated to the later afterglow emission and which is often ex- 
plained as prompt emission originating from higher latitudes 
at the emission source, after the emission at smaller angles 
has faded away. 

In the conte xt of this model, we can test the curvature 
effect a = 2 + (3 (|Kumar fc Panaitescull2000l ). We used the 
values pre- and post-break in the XRT data, showing that by 
this relation, the ob served ai is too shallow and Q2 too steep. 
I Zhang et all (|2006h list several possible causes why this rela- 
tion would not hold. Most important are the effect of using 
the trigger time instead of the time of the last emission peak 
(thereby assigning an incorrect zeropoint for a power law fit), 
and that of superposition of tails of several emission peaks. 
Both causes could be valid here, since the last significant 
emission peaks some 40 seconds after the trigger. However, 
applying such a correction for both effects would effectively 
produce an ever shallower decay, making the discrepancy 
larger. Possible causes for a being too shallow are the f act 
that one is looking at a structured jet l|Rossi et al.ll2002T ) or 
when the observed (X-ray) wavelength regime is still below 
the cooling br eak. In the latter c ase, the decay is given by 
a ~ 1 ± 3/3/2 (jZhang et al-lfeOOq ). Only the very early sec- 
tion of the X-rays agrees with this relation, with a = 0.91 
and f3 = -0.1 (Section [2321 Table[2j). Post-break, however, 
the discrepancy is even larger in this scenario. Potentially, 
both the peak frequency and the cooling break have passed 
the X-ray band at this time. While this could explain the 
relatively long spectral transition, there is no direct evidence 
for this. The discrepancy at later times, for either relation, 
might more likely be attributed to an increasing influence 
of the external shock (which likely dominates the last data) , 
and the last flare, which hamper a good fit to the steeper 
tail of the early light curve. 

Evolution of the early emission in the X-ray to 
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Table 3. The results of our optical follow-up campaign, covering the wavelength range from near infra-red to the near ultraviolet. Other 
data collected through the GCN network provide only upper limits similar to those of the FTN and UVOT, and have therefore not been 
repeated here. 



Start date 


dnys since burst 


exposure time filter 


iiicL^ni tude 


telescope / instrument 


(days UT, July 2005) 


( m id- 1 i me ) 


l LlLLUll LI o I 








16.529991 


0.01954 


400 


V 


> 21.5 


FTN 


16.527721 


0.00306 


30 


R 


> 20.3 


FTN 


16.527721 


0.02239 


540 


R 


> 22.3 


FTN 


16.529911 


0.02236 


520 


B 


> 23.0 


FTN 


16.526169 


0.4315 


5799 


V 


> 21.2 


Swift/VVOT 


16.527859 


0.4042 


5507 


B 


> 22.1 


Swift/VVOT 


16.527697 


0.3988 


5477 


U 


> 21.6 


Swift/VVOT 


16.527534 


0.4416 


5053 


Wl 


> 21.5 


Swift/VVOT 


16.527373 


0.4374 


5486 


Ml 


> 21.8 


Swift/VVOT 


16.528044 


0.4084 


4383 


W2 


> 21.8 


Swift/VVOT 


19.055751 


2.7699 


9480 


'white' 


> 22.5 


Swift/VVOT 


16.579109 


0.05594 


300 


J 


20.76 ± 0.08 


UKIRT/UFTI 


16.581678 


0.05937 


450 


H 


19.31 ± 0.05 


UKIRT/UFTI 


16.564178 


0.04272 


600 


K 


17.74 ±0.15 


UKIRT/UFTI 


16.583727 


0.06055 


300 


K 


18.22 ± 0.15 


UKIRT/UFTI 


17.421337 


0.9275 


3600 


K 


20.5 ± 0.2 


Gemini/NIRI 



Table 4. Summary of the main results obtained in Section [2] We have converted the spectral results from T to 0. 



data section f) a 



X-ray, 103 - 223 seconds 


0.26 ± 0.05 


91 + - 30 

u.yi -o.o9 


X-ray, 263 - 513 seconds 


n sr + 08 
0.8b _ 007 


o 7 o +0.19 
d -'° -0.15 


X-ray, post 1 hour 


i ' ul -0.10 


0.99 ± 0.02 


optical, post 1 hour 


2.94 1 


0.80 t f7 



d = 4.2, therefore no error has been calculated. 



gamma-ray en e rgy r ange has been studied before, eg by 
iFrontera et alj (|200ut ). who looked at the spectral evolution 
for eight BeppoSAX bursts which had been detected by both 
the Gamma-Ray Burst Monitor and one of the Wide-Field 
Cameras. While their findings are similar, an important ad- 
vantage here is that we have an early X-ray light curve 
which, with provisions for some flaring, can be fit with a bro- 
ken power law; the break time of this power law is roughly 
coincident with the moment E p passes through the X-ray 
band. The WFC data generally cover only the very early 
part of the burst, and any underlying (broken) power law 
is not visibile. As a result, the E p evolution could only be 
deduced from the spectral changes for the BeppoSAX data. 



Variable X-ray absorption at early times in the X-ray 
spectra of GRBs has been noted before in other GRBs 
llFrontera et alJ l200d: lAmati et all l200d: lin't Zand et all 
l200ll ; IStarling et alj I2OO0I ). lLazzati fc Pernal l|2002l ) have 
suggested that ionisation strongly modifies the absorption 
properties of the surrounding material. For this effect to 
be noticable, the material should be in a compact region 
surrounding the GRB (< 5 parsec), and have an initial col- 
umn density Ah > 10 21 cm -2 . The change detected here is 
0.8 • 10 21 cm" 2 at redshift 0. For a redshift of z = 2 (see 
below), the column density would be A'h ~ 1.5 • 10 22 cm~ 2 
(at Solar meticallicity) , large enough to fulfill the A'h > 
10 21 cm -2 criterion. 



3.2 The late-time light curve and broad-band 
spectrum 

From the combination of the X-ray temporal and spec- 
tral index around 5ks, we can infer a electron power law 
ind ex of p = 2.0 within the standard fireball model 
(eglMeszaros fc Reedfl997l : ISari et al.lfl999l : IChevalier fell 
1999). This assumes that 1) the afterglow is in the slow- 
cooling regime (i/ m < v c ), 2) i/ c is below the X-rays and 3) 
the jet-break has not manifested itself at this time. Other 
scenarios are easily ruled out by invali dating relations be- 
tween the two indices (see eg Table 1 in lZhang fc Meszarosl 
2004). There is still a degeneracy between the type of cir- 
cumburst environment (a constant density medium or a 
medium with a R~ 2 dependency (7? the distance from the 
centre) such as a stellar wind), which can potentially be re- 
solved by looking at the temporal and spectral indices below 
v c \ in our case the only candidate for that is the optical band, 
which is generally found to be below v c at these times. To 
verify this is indeed the case, we compared the X-ray and 
optical decay indices. 

There is a 1.5cr difference between the afterglow decay 
slope in optical and X-rays, which is only marginally signif- 
icant and does allow v c to be below optical wavelengths. 
However, a cooling break between the optical and X-ray 
waveband makes the picture self-consistent: for an ISM-like, 
the decay index in optical should be 0.25 lower than that 
in X-ray, 0.74, which is certainly compatible with a ob- 
tained from our fit to the 7^-band data. On the other hand, 
for a wind-like medium, the optical decay index would have 
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Figure 5. The X-ray to optical/nIR broad-band spectrum of 
GRB 050716, 0.04 days post burst. The best fit is shown for a 
broken power law connecting the two wavelength ranges, with 
host-galaxy absorption for an SMC-like extinction curve at z = 2. 
The X-ray spectrum has been represented relatively simply by a 
single data point. This has the purpose of enhancing the X 2 e( j 
sensitivity on the nIR data instead of on the X-ray data, since we 
try to calculate the estimated optical/UV host-galaxy extinction. 



to be 0.25 higher than that in X-ray, 1.24, and is therefore 
incompatible with a . 

The non power-law behaviour of the near infra-red 
(nIR) spectrum indicates reddening, which we estimate by 
examining the broad-band spectral behaviour, by combining 
the X-ray, optical and nIR data. For this, we have extrap- 
olated the optical and X-ray data to a common epoch: we 
choose 0.04 days (3.5ks) after trigger, which is central to 
the optical FTN and nIR UKIRT data, and lies just at the 
beginning of the late-time X-ray data. We obtained an unab- 
sorbed X-ray flux from the data between 4ks and and 12ks 
(orbits 2 and 3 combined), which were then extrapolated 
back to 0.04 days using the power law decay index of 0.99. 
The optical and nIR data were corrected for Galactic extinc- 
tion and converted to fluxes as before. We then extrapolated 
to 0.04 days post trigger using our estimate for the optical 
power law decay index. 

With the scaled fluxes, we constructed the broad-band 
spectrum of GRB 050716 at 0.04 days (Figure [5}. We then 
tried to fit the data with models available within the stan- 
dard fireball model, allowing for intrinsic (UV/optical) ex- 
tinction within the host galaxy and a Lyman cutoff result- 
ing from its redshift. No spectroscopic redshift is known, 
however. While in principal one can obtain a photometric 
redshift, there are too many free parameters with the few 
available data to do this accurately. We have therefore ob- 
tained an estimate of the intrinsic extinction for an afterglow 
positioned at various redshifts (1, 2, 3, . . . , 7). We have fur- 
ther constrained the power law slope connecting the X-ray 
and optical/nIR regime to be equal to that obtained by the 
X-ray spectral fit alone, f3 = 1.0. 

A single power law, absorbed by host-galaxy extinction 
or intergalactic Lyman extinction, is not able to reproduce 
the observed broad-band spectrum. A spectral break be- 
tween optical and X-ray is required, as already found from 
the fits to the optical and X-ray light curves. With the as- 
sumption that this break is indeed the cooling break (with 
A/3 = 0.5), we have kept the lower end of this break fixed 



at a power law index of /3 op t = 0.5, and the higher end 
at /3x = 1.0). Reasonable fits are then obtained for red- 
shifts above z = 2, with host-galaxy extinction that follows 
an SMC-like extinction curve (Eb-v = 0.5 - 0.7, Ay = 
1.5 - 2.1, with lower values for higher redshifts. See Figure 
0. This estimate for the host-galaxy is higher than on av- 
erage found for G RB host galaxies, but not unknown (eg 
iLevan et al.ll2006T h However, localised extinction or a line- 
of-sight through an edge-on galaxy could easily explain the 
somewhat larger than average host extinction. 

The reduced \ 2 i s m aJl cases < 1 (for redshifts above 
4, it becomes much less than 1), resulting in a lower limit 
on the redshift of z > 2. A Galactic-type extinction curve 
results in a steeper nIR slope than measured here, and pro- 
vides a worse fit. The same gradual nIR slope also implies 
that host-galaxy extinction should be present, and cannot 
be explained by a Lyman break alone. We further looked 
at the resultant extinction value for small changes (±0.1) 
in our assumed spectral index, which results in changes of 
about 0.1 magnitude in i5s_v(the value for Xrcdi however, 
also increases considerably doing so). Finally, from our best 
fit we find that the v c cannot be located directly below the 
X-ray frequencies, but is restricted to ~ 10 16 Hz. 

An Eb-v of 0.7 with an SMC-like extinction curve 
would imply an iV H of 3.2 ■ 10 22 (|Martin et al.lll989T l in the 
host galaxy. Correcting for an assumed SMC-metallicity of 
1/4 Solar and adjusting for the redshift of z — 2, this turns 
out to be not measurable above the estimated Galactic ex- 
tinction in our data, and indeed, no excess extinction is mea- 
sured past 300 seconds. 

We can obtain a lower limit on the fluence and burst 
energy from the BAT spectrum, which results in 5.9- 10 6 erg 
cm -2 and 5.6T0 52 erg, respectively. Both values are in the 15 
- 150 keV energy range in the observers frame. With a lower 
limit of z > 2, th e relation found by Amati (|Amati et al.l 
l2002tlAmati2006h gives E p > 76 keV in the observers frame, 
while we find E p > 170 keV from the spectral fit (Table [TJ 
after conversion of the cut-off energy to E p ). Note that the 
latter uses a power law with an exponential cut-off, not the 
Band function. Using the previous fit values with a Band 
function, fixing f3 = —2, and extrapolating outside the BAT 
spectral range (0.001 - 10 4 keV), we obtain an estimate of 
the isotropic energy of Ei BO ~ 2.3 ■ 10 53 erg, which is almost 
four times larger than the 15 - 150 keV estimate. Using the 
Amati relation, this would put the peak energy at ~ 150 
keV in the observer frame, assuming z ~ 2. This redshift 
estimate therefore agrees reasonably well with the results 
expected from the Amati relation, and indicates z ~ 2 is a 
likely good estimate for the redshift of this burst. 



4 SUMMARY 

We have presented here a comprehensive analysis of the data 
available on GRB 050716. The most remarkable feature of 
the early emission is the spectral change, clearly visible in 
the X-rays, which can possibly be traced back to the gamma- 
ray emission. The temporal break coincident with the spec- 
tral break points to the peak-flux frequency passing through 
the X-rays at this time. In addition, the X-ray absorption 
before the break appears to be higher than after. While this 
could be an artificial result related to the break, care has 
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been taken to eliminate instrumental effects and the result 
appears to be significant, possibly indicating ionisiation of 
surrounding medium. 

At later times, our multi-wavelength data of 
GRB 050716 shows a fairly standard afterglow: the X- 
ray light curve and spectrum indicate p = 2, with no 
evidence for a break related to the broadening of the 
jet outflow. The optical data are rather limited, partly 
because the optical afterglow is relatively faint: host-galaxy 
extinction with Av = 1.5 - 2.1 can account for the obtained 
upper limits. Though no spectroscopic redshift has been 
obtained for GRB 050716, we obtain a lower limit of z > 2, 
which appears to be in agreement with the Amati relation. 
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